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1  INTRODUCTION  AND  OBJECTIVES 


1 . 1  Introduction 

Materials  having  acceptable  radiation  transmission  properties 
for  use  as  infrared  windows  and  radomes  on  aircraft  generally  have 
poor  resistance  to  fracture  in  particle  impact  situations.  The 
microfractures  that  corm  and  grow  under  bombardment  by  rain,  dust, 
and  ice  particles  scatter  incident  radiation  and  result  in  continu¬ 
ous  degradation  of  the  optical  properties  with  flight  time.  I.onger 
exposures  or  more  severe  impact  environments  can  result  in  interaction 
of  microfractures  with  each  other  and  with  the  radome  surfaces,  causing 
mass  loss  (erosion)  and  gross  mechanical  failure.  Erosion  and  gross 
failure  behavior  at  high  temperature  must  also  be  considered  in 
designing  ceramic  components  of  turbine  engines.  The  Navy  wishes  to 
ascertain  high  temperature  failure  mechanisms,  to  establish  more 
impact-resistant  mierostructural  forms  of  ceramic  materials  and  to 
develop  a  capability  for  predicting  the  fracture  behavior.  In 
support  of  these  needs,  SRI  International  is  performing  research  to 
understand  failure  behavior  of  SijN^  at  turbine  engine  operating 
temperatures,  to  examine  failure  behavior  of  mierostructural ly 
toughened  S^N^  and  to  establish  expressions  and  material  properties 
governing  fracture  development  in  ceramic  window,  radome,  and 
engine  materials  under  hard  particle  impact. 

1 . 2  Obj ect ives 

This  annual  technical  report  documents  the  results  and  progress 
attained  during  the  fourth  year  of  a  research  program  aimed  at 
solving  specific  problems  of  impact  damage  of  ceramics. 

The  three  specific  objectives  in  this  fourth  year  were  to: 

(1)  Determine  the  erosion  behavior  of  Si^N^  at  I400°C 

(gas  turbine  operating  temperature)  and  note  differences 
and  similarities  with  room  temperature  behavior. 
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(2)  Ascertain  ero°ion  performance  of  a  microstructural ly 

toughened  SijN^'ZrCb  ceramic  and  compare  the  results  with 
untoughened  material. 

(3)  Develop  a  predictive  capability  for  calculating  fracture 
damage  in  ceramics  caused  by  particle  impact. 

Completed  accounts  of  the  work  on  objectives  1  and  2  are  presented 
here  as  Sections  2  and  3,  and  have  been  accepted  for  publication  in 
the  Bulletin  of  the  American  Ceramic  Society  and  the  Journal  of 
Materials  Science,  respectively.  The  progress  and  status  of  the 
effort  on  objective  3  are  reported  in  Section  4. 
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2  PARTICLE  IMPACT  DAMAGE  IN  SILICON  NITRIDE  AT  1400°C 


2.1  Background 

The  properties  of  silicon  nitride — high  strength  at  elevated 
temperatures,  low  density,  low  coefficient  of  thermal  expansion,  and 
high  oxidation  resistance — make  it  an  attractive  low-cost  alternative 
to  superalloys  for  stators,  rotors  and  other  high-temperature  engine 
components.  A  major  concern,  however,  with  using  silicon  nitride 
components  in  flying  turbines  is  the  sensitivity  of  their  load-bearing 
capacity  to  flaws  and  cracks.  Aircraft  turbine  engines  ingest  solid 
particles  during  flight,  and  the  resulting  impacts  can  produce  cracks 
that  could  reduce  the  strength  below  that  necessary  to  sustain  design 
loads,  thereby  causing  catastrophic  failure. 

Recent  research  efforts  (Reference  1,  for  example)  have  estab¬ 
lished  the  failure  phenomenology  and  the  strength  degradation  effects--’ 
of  impact-induced  cracks  in  Si^N-i  at  room  temperatures.  However, 
similar  information  is  lacking  at  the  elevated  temperature  representa¬ 
tive  of  gas  turbine  service  conditions.  This  section  describes  initial 
results  of  impact  experiments  on  SijN^  at  1400°C  and  compares  the 
damage  observations  with  those  from  room  temperature  experiments. 

2 . 2  Experimental 

Octagonal  specimens  roughly  25  mm  in  diameter  were  cut  from  a 
b . 55— mm— thick  plate  of  fully  dense  NC132  hot  pressed  silicon  nitride 
(HP  Si^N^)  purchased  from  the  Norton  Company,  Worchester,  Mass.  The 
microstructure  consisted  of  '3-SijN^  grains,  most  of  which  were  equi- 
axed  and  less  than  1  vim  in  size,  although,  some  were  up  to  4  urn  long 
with  aspect  ratios  up  to  10.  No  significant  grain  boundary  or  in¬ 
clusion  phases  were  observed,  but  traces  of  tungsten  carbide  and 
tungsten  disilicide  were  detected  by  x-ray  diffraction. 

The  impact  side  of  the  specimens  was  polished  and  the  specimens 
were  placed  into  a  central  vertical  refractory  tube  of  a  furnace  onto 
a  12-mm-thick  silicon  nitride  slab  supported  by  an  inner  tube.  The 
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particle  impact  facilitv  is  shown  in  Figure  1.  To  reduce  possible 
effects  of  stress  wave  reflections  from  the  rear  specimen  surface, 
the  mating  surfaces  of  specimen  and  slab  were  also  ground  flat  and 
polished.  The  temperature  of  the  specimens  was  increased  to  1  iBO  - 
20°C  in  about  3  hours  and  maintained  at  this  level  for  the  particle 
impact  experiments.  A  platinum-rhenium  thermocouple  monitored  tin- 
temperature  at  the  specimen  location.  While*  they  were  at  elevated 
temperature,  the  specimens  were  bathed  in  a  stream  of  nitrogen  gas 
to  reduce  surface  oxidation. 

Single  spheres  of  tungsten  carbide  or  steel  2.4  mm  in  diameter 
were  accelerated  in  a  small  gun  barrel  by  the  sudden  release  of  com¬ 
pressed  air.  A  thin  disk  located  between  the  gas  reservoir  and  the 
gun  barrel  sustained  the  compressed  air  until  a  pressure  was  reached 
that  caused  the  disk  to  rupture.  Rupture  pressure,  and  hence  particle 
velocity,  was  varied  by  using  disks  of  different  materials  and  thick¬ 
nesses.  Velocities  below  85  m/s  were  achieved  by  omitting  the  rupture 
disk,  pressurizing  the  compressed  air  reservoir  to  a  desired  low 
pressure,  and  opening  the  valve  by  hand.  Impact  velocities  were 
measured  by  the  photomultiplier  arrangement.  Light  reflected  by  the 
particle  into  a  photomultiplier  at  three  fixed  stations  indicated 
particle  position  at  three  times  during  its  flight  to  the  target. 

2 . 3  Resu 1 ts 

Room  temperature  experiments  were  carried  out^  at  velocities  ranging 
from  7  to  230  m/s.  No  damage  was  produced  by  impacts  with  2. 4 -nun-diameter 
tungsten  carbide  spheres  at  velocities  less  than  15  m/s.  A  partial  ring 
crack  appeared  at  15  m/s.  Velocities  between  15  and  30  m/s  produced  ring 
cracks  of  increasing  number  and  size,  and  velocities  between  30  and  90  m/s 
produced  ring  cracks,  cone  cracks,  and  a  plastic  impression.  Above  90  m/s, 
radial  cracking  as  well  as  fragmentation  of  the  tungsten  carbide  sphere 
occurred . 

Impact  experiments  at  1400°C  were  performed  at  velocities  of  38,  51, 
85,  121,  156,  173,  and  197  m/s  with  tungsten  carbide  spheres.  Radial 
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tracks  and  a  well-defined  plastic  imprest'-  .  were  produced  in  each  exper¬ 
iment.  As  shown  in  Figure  _  ,  the  size  of  the  impression  and  intensity 
of  tin  .  .njial  cracks  increased  with  increasing  velocity.  At  velocities 
greater  than  150  m/s,  the  radial  c racks  extended  to  the  sides  and  rear 
surt  aces  ot  the  specimens,  often  fragment ing  them.  Figure  1  compares 
the  damage  at  20l C  and  1400°C  produced  by  tungsten  carbide  spheres  im¬ 
pacting  at  about  50  m/s.  The  type  of  damage  at  14001 C  (radial  cracks) 
is  more  deleterious  than  the  ring-cone  cracking  that  occurs  at  20°C  be¬ 
cause  radial  cracks  tend  to  be  longer  and  deeper.  Moreover,  the  damaged 
surface  area  is  substantially  larger  at  1400°C.  The  relative  dynamic 
hardness  at  1400°C,  computed  from  the  plastic  impression  diameters. 

Figure  2,  using  the  Brinell  equation,  decreased  mono ton i cal  1 v  from  0.8 
at  50  m/s  to  about  0.6  at  170  m/s. 

Cone  cracks  occurred  along  with  radial  and  lateral  cracks  in  Si^N, 
at  1400LC,  as  evidenced  by  the  cross  section  shown  in  Figure  4.  This 
fracture  surface,  produced  by  extension  of  radial  cracks,  reveals  the 
site  of  the  cone-shaped  fragment  beneath  the  impact  site,  the  lateral 
cracks  that  initiate  internally  near  the  boundary  of  the  plastic  zone, 
the  profile  of  the  radial  cracks,  and  the  mic  rocrac’-.ed  zone  direct  lv 
beneath  the  plastic  impression. 

By  contrast,  the  fracture  damage  produced  by  2 . 4-mm-d iameter  steel 
spheres  consisted  of  annuli  of  segmental  ring  cracks  both  at  room  tem¬ 
perature  and  at  1400°C,  even  at  velocities  of  230  m/s,*  Figure  5.  The 
threshold  velocity  for  fracture  at  20°C  by  steel  spheres  was  about  40  m/s, 
or  nearly  three  times  as  high  as  the  threshold  velocity  for  tungsten 
carbide  spheres. 

2.4  Discuss i on 

These  results  show  that  the  impact-induced  fracture  pattern  in 
Si.^N,  it  1400 can  he,  but  is  not  necessarily,  different  than  that  at 
20°C.  The  beh  avior  seems  to  depend  on  the  hardness  of  the  particle 
relative  to  the  target,  and  can  be  explained  by  a  difference  in  the 
near-crater  stress  distribution  when  elastic  or  plastic  deformation 
occurs  in  the  target. 
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FIGURE  2  PLASTIC  IMPRESSION  DIAMETER  VERSUS  IMPACT  VELOCITY 
FOR  HP  S13N4  AT  ROOM  TEMPERATURE  AND  1400  C 

(2.4-mm-dia.  tungsten  carbide  spheres) 
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FIGURE  4  SIDE  VIEWS  OF  SPECIMEN  FRACTURED  BY  AN  IMPACTING  TUNGSTEN  CARBIDE  SPHERE 


have  explained  tin-  various  crack  i  up, 
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Evans,  Lawn,  and  Wilsliaw  ’ 
patterns  observed  in  ceramics  around  indenters  and  impacting  particles 
on  tile  basis  ol  dittercnces  in  stress  states  caused  by  the  absett  c  or 
tile  occurrence  of  plastic  flow.  When  the  target  deforms  *■  1  ast  i  ea  I  I  v , 
tlie  maximum  tensile  stresses  are  directed  radial Iv,  resulting  in  an 
annulus  of  segmented  ring  cracks  just  outside  the  c i re  I c  of  cont  it  t. 

Higher  indentation  pressures  or  impact  velocities  cause  these  shallow 
ring  cracks,  which  grow  initially  nearly  normal  to  the  , out  a.  t  surtnee, 
to  extend  downward  and  away  ‘rom  tlie  zone  of  compression  under  the  ,  on- 
tact  site  to  I  orni  cone-shaped  cracks  with  sides  rough  lv  40  to  tin  imp.,,  t 
surface  or  22*  to  the  indentation  surface. 

When  the  material  under  the  contact  area  deforms  plastics!  !v,  radial 
stresses  relax  and  the  stresses  in  the  tangential  direction  he,  o:m-  let:  i).  , 
often  resulting  in  the  appearance  of  radial  ,  racks  extending  outward  ;  t 
the  periphery  of  the  plastic  impression  in  a  spokelike  array.  Tcnsil* 
stresses  also  arise  in  the  specimen  interior  normal  to  the  impacted  stir!  «,  <  , 
particularly  during  unloading,  and  may  produce  lateral  <  racks  parall,  I  t  ■ 
the  specimen  surface. 

These  considerations  suggest  that  when  the  Si^N^  is  at  r*00‘V,  tin 
relative  hardnesses  of  Si.,N,  and  tungsten  carbide  may  he  reversed,  a  1  1  >>\.  i  u .• 

)  *4 

the  tungsten  carbide  sphere  to  remain  elastic  and  transmit  loads  to  tin 
Si  N,  target  that  exceed  the  high  temperature  vicld  strength  ot  Si  N,. 

J  *4  ‘  >  4 

In  tlie  vast'  of  steel  spheres,  if  the  hardness  of  Si.^N,  at  1*00  ('.  remain- 
greater  than  the  room  temperature  hardness  of  the  steel,  then  when  impn  t 

occurs,  the  sphere  would  deform  plastically,  preventing  loads  in  excess 

of  the  Si.,N.  vield  strength  and  resulting  in  elastic  fracture. 

J  4  ‘ 

These  inferences  are  in  accord  with  relative  dynamic  hardnesses 
estimated  t rom  static  room  temperature  hardness  values  of  Ih(lt),  1470, 
and  820  kg/mm^  measured  tor  the  silicon  nitride  plates,  the  tungsten 
carbide  spheres,  and  the  steel  spheres,  respectively.  The  impression 
diameter  data  in  figure  2  show  that  the  hardness  of  Si-jN,  at  I  *00 "C 
is  about  807'  of  the  room  temperature  hardness,  or  about  1100  kg/mnr. 
Therefore,  assuming  similar  rate  effects  for  Si  and  WC,  the  WC 
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spheres  should  i  .nisi'  clast  io  fracture  behavior  in  Si  jN^  at  mm  tempi- r- 
aturc,  but  c  I  as  t  i  c  /  p  I  as  t  i  c  tract  tire  behavior  at  l'tOO"C,  as  was  observed, 
f ti r t hormore ,  since  the  hardness  o!  the  steel  spheres  is  tar  below  the 
Si  hardness  even  at  ItOU'C,  steel  spheres  should  produce  on  1  v 
elastic  I  racture  ot  the  Si  j\(,  a  1  so  as  observed .  The  slight  differences 
in  hardness  ot  Si  jN  (  and  WC  (lht)i)  versus  I  a  70  kp/min-')  will  likelv  be 
enhanced  in  these  experiments  hv  the  diitcrcnces  in  geometry  for  the 
two  materials.  (leometries  that  constrain  the  plastically  deformed  zone 
mo  re  h  i  yh  I  v  (such  as  the  Si  jN  plate,  wh  i  ■  h  surrounds  the  p  1  ,,s  t  i  < 
impression  hv  a  halt  spaa  ot  hulk  elastic  material)  will  require 
higher  loads  tor  a  pi  veil  indentation  than  less  constrained  pcometries5 
(such  as  thi'  smu  I  i  impacting  sphere). 

Studn.an  and  field  observed  both  radial  and  rinp  crackinp  in  a 
quenched  steel  arcuiti!  indentations  produ.  ed  hv  tun  vs  ten  carbide  spheres 
and  also  explained  tile  behavior  in  terms  ot  a  '  i  ansi  t  ion  trot'  clast:, 
to  c  1  asl  ie/p  1  as  t  i  c  behavior.1’ 

-  .  ">  (lone  1  us  i  ons  and  Ke>  ommeiidal  ions 

These  results  point  to  puidt  lines  tor  predicting  the  type  and  extent 

of  damape  in  Si  N,  turbine  components  to  be  expeited  t  run  riven  dibris 
\  -* 

particles.  It  indeed  the  felativi  liardtiesses  ot  parti.  1«  and  .  ompniunt 

govern  the  tvpe  ot  crackinp  in  tin  component,  a  know!. dpt  ot  tin-  component 

hardness  at  tin  t  emper.it  lire  ot  .ervii  e  and  the  hardness  ot  tin  ant  i'  ipated 
debris  particles  wouid  indicate  tiie  I rai  turc  peomet rv  and  allow  the  pro- 
pensitv  for  streiipth  deprad.it  i  on  and  erosion  to  be  asst  ss>  d .  (r.nks  that 
tend  to  penetrate  the  lompotietit  (radials,  medians  and  tones!  1 1  ml  to 
tlepratle  its  stretipth,  whereas  cracks  that  t  t  nil  to  turn  ha.  k  to  tin  -uir  i  a.  a 
(laterals  and  tones.)  tend  to  promott  mass  loss  and  ir.-ion.  A  "tor.  qtian- 
t  i  t  a  t  i  ve  predictive  c  a  pah  i  I  i  t  v  must  aw  a  i  t  add  i  t  i  on  a  1  .  xp.  I  i  men  t  ■;  tii.M  show 
the  extent  of  dnmape  and  the  result  inp  st  n  net  It  tit  p  r.uia  i  i  on  and  erosion 

lit  ha  v  io  r  at.  a  !  into  t  i  on  ot  relative  t  a  r  p.e  i  and  p.t  r  t  i .  1  <  ha  rtliu  '■'•!  s ,  part  it  1 . 

velocity,  sire,  shapt  ,  anti  the  number  ol  impacls. 


3  IMPROVLI)  IMPACT  I'RACTl'KK  KKS 1 STANC1.  IN 
OX  1  DAT  1  ON-Torcill'NKl)  Si-N/.-dt)  VOL  Zrl>  > 


3.  1  Background 

An  oxidation  technique  to  induce  compress  i  vo  st  ri'sscs  in  the 
ni'.ii — surface  regions  of  Si  jN/4-ZrO}  spi-r  imens  lias  been  repot  tod  liv 

-J 

Lange.  '  The  material  is  made  bv  thoroughly  mixing  si  I  iron  nitride 
powder  with  about  20  vole  zirconium  oxide  and  A  vo  I  ‘  A1>0.  powders, 
then  hot -press  i  np  to  achieve  .i  fullv  dense  plate.  Ihc  plate  is  then 
given  a  subsequent  anneal  at  700°C  for  3  hrs.  in  air.  The  zirconium 
oxynitride  serond-phase  within  several  grain  diameters  ot  the 
surfaces  is  oxidized  to  monoclinic  Z  rO  j  with  an  accompany i ng  increase 
in  molar  volume  of  about  •»->?.  The  expansion  of  the  lattice  in  tin- 
surface  regions  gives  rise  to  substantial  compressive  stresses  on  tin 
surface  (and  significant  tensile  stresses  in  the  interior).  Tin- 
oxidation  increases  the  apparent  surface  toughness  and  flexural 
strength  by  about  25%.?  These  results  encouraged  us  to  investigate 
the  response  of  this  material  to  particle  impact. 

3.2  Materials,  experiments,  and  Result  s 

Lange  provided  two  specimens  each  of  oxidized  and  unoxidized 
Si  jM/,-20  vo  1 7.  ZrO)  as  halves  of  broken  4-point-bend  specimens.  Nominal 
dimensions  were  1  x  h  x  18  mm.  The  surface  finish  ol  the  unoxidized 
specimens  was  improved  by  polishing  to  make  it  comparable  to  tin- 
surface  finish  of  the  oxidized  specimens. 

Using  1 . 2-mm-d  i  atiii-t  i-r  tungsten  carbide  spheres  and  a  compressed 
air  gun , e  we  conducted  a  series  of  impact  experiments  at  various 
velocities  on  oxidized  (surface-strengthened)  and  unoxidized 
S  i  )N/(-Zrt)  i  to  determine  the  damage  phenomenology,  the  threshold 
velocities  tor  ring  and  radial  cracks,  and  the  rates  ol  impact  damage 
deve 1 opment . 

At  the  lowest  velocities,  no  damage  could  In-  detected  at 
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magnifications  up  to  400X.  Slightly  higher  velocities  produced  a 
plastic  impression,  and  at  higher  velocities,  ring  cracking  began, 
followed  by  radial  cracking  at  still  higher  velocities. 

The  size  of  the  plastic  impression  increased  mono! on i ca 1 1 v  with 
increasing  impact  velocity.  A  prof i lometer ,  used  to  measure  the 
shapes  of  the  impression,  verified  that  plastic  f 1 ow  at  the  impact 
sites  occurred  at  velocities  well  below  those  necessary  for  ring 

"k 

crack  formation.  At  low  impact  velocities,  larger  impressions  wire 
sustained  by  oxidized  material.  Figure  6,  showing  thaL  the  oxidation 
treatment  decreased  the  dynamic  hardness  of  the  surface  layers.  At 
higher  velocities  (>80  m/s)  impression  diameters  of  oxidized  and 
unoxidized  material  are  nearly  equivalent,  indicating  that  the  soften¬ 
ing  effect  is  confined  to  the  surface  layers. 

A  similar  effect  on  quasi-static  hardness  was  found  in  a  series 
of  diamond  pyramid  hardness  indents  at  loads  ranging  from  0.1  to  20  kg. 
The  average  hardness  of  the  oxidized  material  was  1200  kg/mnr, 
compared  with  a  value  of  1600  kg/mm^  for  unoxidized  specimens.  Again, 
a  more  pronounced  softening  in  the  surface  layers  of  the  oxidized 
specimens  relative  to  the  bulk  was  indicated  by  the  hardness  results 
(1000  kg/mnr)  at  0.1  and  0.2  kg  indent er  loads. 

The  threshold  velocity  for  fracture  damage  in  the  oxidized 
material  (~45  m/s)  was  substantially  higher  than  for  the  unoxidized 
material  (~1 7  m/s),  and  the  amount  of  damagi  produced  by  a  given 
higher  velocity  impact  was  significantly  less.  At  an  impact  velo»  itv 
of  45  m/s  a  tort  ion  of  a  single  ring  crack  was  produced  in  the  oxidized 


*Tbis  behavior  is  opposite  that  observed  bv  A.  G.  Ivans,  I'niversity 
ol  California,  Berkeley  (private  commun i  ■  at i on ) ,  and  i  ons i d<  a  a t  i oils  el 
tile  stress  distribution  around  a  plastic  contact  show  that  radial 
tracks  indeed  should  be  expected  because  el  large  ,  i  iv  unit  e  r<  n  t  i  a  1 
tensile  stresses.  blast  ic  contact  prodtn  es  large  stresses  in  the 
radial  direction  and  hence  encourages  ring  cracking.  In  the  elast  i. - 
plastic  regime,  however,  the  stress  distribution  nav  t  aver  ring,  ot 
radial  cracks,  and  thereby  account  lor  tin  d i sc  repair  \  in  observed 
f  rue t  ur i ng. 
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FIGURE  6  PLASTIC  IMPRESSION  DIAMETER  VERSUS  IMPACT  VELOCITY 
FOR  OXIDIZED  AND  UNUXIDIZED  Si  ,N4  20  2rO, 


material,  whereas  several  wel 1 -developed  ring  cracks  were  produced  in 
the  unoxidized  specimens.  Ring  cracks  became  more  numerous  as 
velocity  increased.  Figure  7  shows  the  damage  produced  in  the  two 
materials  by  145  m/s  impact.*  An  annulus  of  cracked  material  deve¬ 
loped  about  the  center  of  the  impact  site.  The  annular  area  was 
always  greater  in  the  unoxidized  material.  At  a  velocity  of  180  m/s, 
radial  cracks  appeared  in  the  unoxidized  material.  No  radial  cracks 
were  observed,  however,  in  oxidized  material  at  velocities  up  to  195 
m/s.  Figure  8  shows  that  the  inner  radius  of  the  cracked  zone  is 
larger  for  the  oxidized  material,  and  the  outer  radius  is  smaller. 
Thus,  for  any  given  impact  velocity,  the  damaged  area  in  oxidized 
material  is  significantly  smaller  than  in  the  untreated  material. 

Several  impact  sites  in  oxidized  and  unoxidized  material  were 
sectioned,  and  the  section  surfaces  were  polished  and  examined  with 
a  microscope  to  determine  the  cracking  patterns  in  the  specimen 
interior.  Both  materials  showed  the  Hertzian  cone-shaped  cracks. 

3. 3  Discussion 

The  observed  enhancement  of  impact  damage  resistance  could  be 
attributable  to  oxidation-induced  surface  compressive  stresses  or  to 
oxidat ion- induced  softening  or  to  both.  However,  the  important 
aspect  of  the  impact  damage  resistance  is  in  the  fracture  propensity 
of  the  two  materials.  Ring  cracks  and  radial  cracks  were  initiated 
at  lower  Levels  and  developed  at  more  rapid  rates  in  the  unoxidized 
material,  which  suggests  that  the  oxidation-induced  surface 
compressive  stresses  inhibit  and  suppress  fracture. 

Although  the  experimental  results  indicate  superiority  of  the 
oxidized  material  in  resisting  damage  from  an  impacting  particle,  the 
used  of  oxidized  material  in  high-temperature  oxidizing  environments. 


*The  halo  that  is  evident  in  Figure  7  is  shallow  surface  damage 
produced  when  the  tungsten  carbide  sphere  fragments  and  impacts  the 
surface.  Tungsten  carbide  spheres  break  at  velocities  ot  about 
100  m/s  and  greater  when  impacted  against  these  specimens. 
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FIGURE  8  GROWTH  OF  DAMAGE  ZONE  WITH  IMPACT  VELOCITY 
FOR  OXIDIZED  AND  UNOXIDIZED  Si3N4-20%  vol  %  Zr02 
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such  us  those  encountered  by  turbine  blades,  appears  of  doubtful  benefit 
because  previous  impact  tests  indicate  that  over-oxidation  increases 
erosion  rates.' 

We  note  that  the  material  used  here  was  probably  not  optimal; 
other  compositions  and  oxidation  treatments  could  be  expected  to 
exhibit  greater  dynamic  fracture  resistance.  Additional  experiments 
on  other  oxidation-toughened  materials  and  with  particles  of  other 
sizes,  shapes,  types,  angles  of  incidence,  and  velocities  are  suggested 
to  confirm  and  further  investigate  these  promising  initial  results. 
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PROGRESS  IN'  DEVELOPING  A  PREDICTIVE  CAPABILITY 
FOR  HARD  PARTICLE  IMPACT  DAMAGE  IN  CERAMICS 


Development  of  a  quantitative  predictive  capability  for  particle 
impact  damage  in  ceramics  requires  a  detailed  knowledge  of  the  stress 
field  histories  experienced  by  material  in  the  virinitv  of  the  impart 
site  and  mathematical  expressions  describing  the  development  of 
plastic  deformation  and  cracking.  This  section  describes  the 
procedures  and  results  of  a  method  we  have  developed  for  computing 
stress  field  histories  and  presents  the  rudiments  of  a  computational 
fracture  model. 


4.1  Approach 

Our  approach  in  developing  a  predictive  capability  for  particle 
impact  damage  in  ceramics  consists  of  the  following  tasks: 

(1)  Establish  an  experiment  that  can  be  analyzed  to  allow 
identification  and  evaluation  of  the  material  properties 
governing  resistance  of  semibrittle  materials  to  fracture 
damage  by  particle  impact.  Experiments  having  the  high¬ 
est  probability  of  success  appeared  to  be  those  closely 
simulating  particle  impact  conditions,  but  having  a  simple 
impact  geometry  to  aid  in  interpretation  and  analysis. 

Thus,  we  performed  experiments  involving  normal  (90°) 
impact  of  a  single,  elastic  sphere  against  a  flat  specimen 
sur face . 

(2)  Obtain  an  understanding  of  how  fracture  damage  develops 
about  an  impact  site.  We  studied  damage  phenomenology 
by  performing  experiments  at  various  velocities  so  as  to 
produce  various  levels  of  damage  ranging  from  below 
incipient  to  surface  fragmentation.  By  examining  the 
impact  sites  and  cross  sections  taken  through  the  sites 
with  a  microscope  we  deduced  how  fracture  damage  evolves. 

(3)  Describe  quantitatively  the  deformation  and  fracture 
damage  produced  by  particle  impact.  Definition  of  impact 
erosion  properties  requires  a  knowledge  of  the  relationships 
between  loading  parameters  and  material  failure  response, 
and  therefore  a  quantitative  description  of  deformation 
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and  f  rarluri'  is  requ i  red .  We  i  mmtud  anil  measured  tilt- 
radial,  ring,  and  lateral  travks  at  particle  impart  sites 
to  obtain  trank  size  distributions  tor  various  impart 
vr 1  or i t its. 

(4)  Computt  stress  field  histories  produced  in  materials  impart¬ 
ed  by  a  solid  particle.  The  stresses  and  strains  and  their 
variation  with  time  and  location  in  the  target  must  In- 
known  in  carder  to  analyze  the  observed  flow  and  fracture 
damage.  Further,  a  procedure  for  computing  stress  field 
histories  under  particle  impact  is  necessary  as  the  basis 
for  the  predictive  capability. 

(5)  Construct  a  mathematical  model  of  the  fracture  processes 

by  making  use  of  fracture  mechanics  concepts  and  correlations 
of  fracture  damage  with  computed  stress  histories.  Write  the 
model  as  a  subroutine  for  the  computational  code  used  in 
the  previous  task  to  calculate  the  stress  field  history. 

(6)  Verify  the  fracture  model  by  using  it  in  the  stress  field 
history  code  to  simulate  selected  particle  impact  experiments. 
Compare  computed  and  observed  fracture  damage,  and  modify 

the  model  until  good  agreement  is  obtained.  Extract  f^om 
the  model  the  expressions  and  properties  describing  impact 
damage  development. 


Tasks  1  through  3  have  been  completed  and  the  procedures  and 
results  are  documented  in  References  1  and  8.  Tungsten  carbide 
spheres  were  accelerated  in  a  pneumatic  gun  to  velocities  ranging 
from  10  to  350  m/s  and  made  to  perpendicularly  impact  polished 
surfaces  of  chemical-vapor-deposited  zinc  sulfide  and  hot  pressed 
Si^N^.  The  individual  impacts  produced  various  levels  of  damage 
in  the  form  of  plastic  impressions  and  populations  of  several  types 
of  cracks.  The  fracture  damage  was  quantified  by  counting  and 
measuring  individual  cracks. 

Damage  threshold  conditions  and  characteristic  rates  for  crack 
nucleation  and  growth  were  deduced  from  the  data.  The  shapes  of  the 
crack  size  distribution  curves  suggest  that  the  size  distribution  of 
inherent  flaws  is  exponential,  consistent  with  the  results  of 
probability  theory,  and  show  that  the  growth  law  for  radial  cracks 
is  similar  in  form  to  expressions  derived  from  elastic  energy  balances. 
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The  next  step  (Task  4)  in  developing  the  predictive  capability  is 
to  compute  the  stress  field  histories  experienced  by  ceramic  plates 
impacted  by  a  sphere,  so  that  the  quantitative  fracture  data  can  be 
correlated  with  the  causative  stress  pulses  (Task  5).  During  this 
reporting  period  the  stress  field  history  computation.  Task  4,  was 
completed  and  fracture  model  development.  Task  5,  was  begun.  This 
section  presents  the  results  of  this  work. 

4.2  Calculation  of  the _ Stress  Field  Jlistorv_  in  a  Plate  Impacted  by 

a  Sphere 

4.2.1  Background 

Computation  of  the  stress  field  history  produced  in  a 
target  material  by  an  impacting  particle  is  complicated  and 
quite  formidable,  involving  elastic  and  plastic  deformation, 
dynamic  effects,  and  the  stress  relaxation  effects  of  cracks 
as  they  deveLop.  Analytical  solutions  exist  only  for  the  most 
simple  (and  least  interesting)  of  situations — namely,  the 
static  indentation  elastic  case  in  which  no  fracture  orcurs.^’^ 
Even  in  this  case,  however,  there  appears  to  be  no  simple  single 
parameter  that  characterizes  the  stress  field;  that  is,  there 
is  no  parameter  analogous  to  the  stress  intensity  factor  that 
characterizes  the  magnitude  of  the  stress  field  at  the  tip  of 
a  crack.*  Such  a  parameter  would  clearly  be  very  useful  in 
correlating  stress  field  history  with  fracture  damage  and, 
hence,  was  sought  during  this  present  research  year. 

We  considered  therefore  numerical  methods  and  examined 
the  various  available  advanced  computer  codes  to  determine  (1) 
if  the  stress  field  history  in  a  target  impacted  by  a  hard 

*The  nonlinearity  arises  because  the  contact  area  increases  non¬ 
linear  ly  with  pressure. 
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sphere  had  already  been  computed,  and  if  not  (2)  whether  such  a 
computation  appeared  possible  with  an  existing  code.  We  found 
that  the  water  drop  impingement  problem  had  been  simulated  by 
Rosenblatt,  Ito,  *1*12  ancj  others  and  that  computations  of 
impact  by  penetrators  with  noses  of  various  shapes  (including 
hemispherical)  had  also  been  completed.  Rosenblatt  and 
Egguml*4  appear  to  have  performed  the  only  computation  of  hard 
sphere  impact.  These  workers  simulated  the  impact  of  a  400- 
pm-dia.  WC  sphere  onto  a  flat  ZnS  surface  at  850  m/s.  They 
assumed  the  sphere  to  be  rigid,  the  ZnS  to  be  elastic-perfect ly 
plastic  and  to  sustain  no  fracture,  and  the  sphere-target  inter¬ 
face  to  be  frictionless.  Under  these  assumptions  they  obtained 
the  contact  pressure  and  contact  radius  as  a  function  of  time, 
the  growtli  of  the  plastic  zone  beneath  the  contact  site,  and 
the  magnitude  and  directions  of  the  tensile  stresses  outside 
the  plastic  zone.  Although  their  results  will  be  modified  by 
the  fracture  activity  that  occurs  in  both  the  target  and  the 
sphere,  the  results  are  useful  in  interpreting,  at  least  quali¬ 
tatively,  the  observed  fracture  morphology  in  ZnS.-*'’  We  therefore 
proceeded  on  a  similar  path,  using  a  2-D  Lagrangian  wave  pro¬ 
pagation  code  to  calculate  the  deformation  and  stress  histories 
in  our  targets.  Development  of  a  reliable  predictive  capability, 
however,  requires  that  the  stress  relaxation  caused  by  the 
developing  damage  be  accounted  for.  We  are  therefore  also 
attempting  to  obtain  a  treatment  for  stress  relaxation  by 
appropriately  modeling  the  fracture  process,  Section  4.3. 

4.2.2  Comput ational  Pro c e  du re 

The  particle  impact  experiments  were  simulated  computation¬ 
ally  using  TROTT,  a  two-dimensional,  Lagrangian,  finite  difference, 
wave  propagation  code  written  at  SRI  Intc  ional.^1  The 
calculational  procedure  is  the  standard  leapfrog  method  of 
Von  Neumann  and  Richtmyer,^  using  artificial  viscosity  to 
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smooth  shock  fronts  over  several  computational  cells,  which 
may  be  quadrilateral  or  triangular  and  treated  as  finite 
elements  for  mass  and  momentum  calculations.  The  code  can 
accommodate  complicated  material  models  because  of  its  large 
storage  capability  per  cell,  but  the  code  is  simple  and  easily 
modified.  Because  of  its  simplicity,  TROTT  runs  about  three 
times  faster  than  standard  2-D  codes. 

In  preparation  for  the  simulation,  a  slide  line  was 
inserted  at  the  particle-target  interface  to  allow  the  spherical 
impactor  to  start  out  not  in  contact  with  the  target,  then  to 
penetrate  an  appropriate  distance  into  the  target  with  sliding 
occurring,  and  finally  to  separate  from  the  target  as  the 
sphere  rebounds.  The  surface  of  contact  was  modeled  by  a 
standard  master/slave  slide  line  (with  both  force  and  displace¬ 
ment  condition)  across  the  sphere- target  interface,  before 
impact  and  after  rebound,  when  the  two  surfaces  are  not  in 
contact,  the  forces  and  displacements  in  the  sphere  and  targi t 
are  calculated  independently.  During  contact,  however,  the 
force  on  both  sides  of  the  contact  surface  govern  the  motion  of 
the  master  side  of  the  slide  line,  and  the  displacement  of  the 
master  side  of  the  slide  line  governs  the  normal  displacement 
of  the  slave  side.  Hence,  the  slave  surface  slides  on  the 
master  surface. 

The  computational  grids  of  the  target  and  spheres  were 
constructed  to  conform  to  a  Cartesian  coordinate  system  and  to 
provide  higher  resolution  in  the  area  of  contact  where  stress 
and  strain  gradients  are  steep.  The  target  cells  were  nearly 
square  and  similar  in  size  to  neighboring  cells.  The  largest 
cell  dimension  did  not  exceed  the  smallest  dimension  by  more 
than  10%.  Likewise,  the  areas  of  neighboring  cells  did  not  vary 
by  more  than  10%.  The  sphere  was  zoned  by  laying  out  a  grid 
of  horizontal  and  vertical  lines  that  had  the  necessary  rosoluti 
and  cell  size  variation,  rotating  this  network  90°  and  mapping 


the  network  into  a  circle  by  means  of  the  Joukowski  trans¬ 
formation.^®  This  transformation  7/  =  7.  +  b^/Z  in  the  complex 
plane  maps  circles  into  lines  and  is  well  known  in  the  field 
of  fluid  mechanics.  The  resulting  grids  are  shown  in  Figure  9. 

We  chose  to  simulate  the  impact  of  an  800-um-dia.  tungsten 
carbide  sphere  with  a  plate  of  CVD  ZnS  at  velocities  of  5,  30, 
and  60  m/s.  Cracking  around  the  impact  site  and  the  concurrent 
stress  relaxation  were  not  treated  in  the  initial  simulations. 
These  important  effects  will  be  taken  into  account  in  subsequent 
calculations  by  the  dynamic  fracture  model  being  developed  in 
Task  5  of  this  research  program.  The  tungsten  carbide  sphere 
was  treated  as  elastic,  but  plastic  flow  of  the  target  was 
allowed  to  occur  according  to  a  stress-strain  relation  deduced 
from  hardness  tests.  Using  the  method  of  Evans  and  Wilshat. 
we  characterized  the  quasistatic  work  hardening  behavior  of  ZnS 
by  plotting  in  log-log  space  the  indentation  pressure  versus 
the  normalized  plastic  impression  diameter  d/D,  where  d  is  the 
impression  diameter  and  D  is  the  diameter  of  the  WC  sphere. 

Figure  10.  The  slope  was  taken  as  an  approximate  value  of  the 
work  hardening  exponent,  n  =  0.3.  The  material  properties  used 
in  the  simulation  are  listed  in  Table  I. 

4.2.3  Results 

The  computation  simulated  qualitatively  all  the  important 
aspects  of  a  particle  impact  event.  The  sphere  approached  the 
target  at  the  prescribed  velocity,  made  contact,  and  began  to 
penetrate.  The  stresses  in  the  target  increased  monoton ica I  I v 
and  when  the  yield  strength  was  reached,  material  in  the 
immediate  vicinity  of  the  impact  site  began  to  flow  plastically. 

As  the  elastic  sphere  continued  to  penetrate,  the  elastically 
deformed  zone  expanded. 

Sphere  velocity  decreased  monoton ica 1 1 v  with  increasing 
penetration,  reached  zero,  then  became  negative  (i.o.,  rebounded). 
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For  clarity,  the  qnd  shown  m  (a)  is  twee  <is  coarse  as  actual  qntf  istmw  »»  fh'l 

All  dimensions  m  millimeters,  MA  4928- 10H 


FIGURE  9  COMPUTATIONAL  GRIDS  USED  IN  SIMULATING  A  SPHERE  IMPACTING  A  PLATE 
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FIGURE  10  STRESS-STRAIN  CURVE  FOR  ZnS  INDENTED 
WITH  AN  800-um-dia.  WC  SPHERE 


at  i'i‘  li-r.it  ing  in  tin.'  direction  opposite  tin'  impact  directum  and 
finally  reaching  a  s t eadv-s t at e  rebound  vc  1  oi  ity  as  it  lost 
cunt. Ht  with  the  target.  Crater  depth  .and  diameter  relaxed 
somewhat  from  the  maximum  values  attained  at  deepest  sphere 
penetration,  and  a  residual  stress  field  remained  about  the 
i mpae  t  site. 

To  verilv  the  quant i t at i Ve  tidelitv  of  the  eomputation, 
we  compared,  wherever  possible,  computed  results  with  measurements 
we  made  in  the  experiments.  I  ah  I e  II  compares  computed  and 
observed  sphere  rebound  velocities,  crater  diameters,  and 
crater  depths  at  three  conditions  o!  impact  velocity. 

The  cornput a t i onn 1  results  agree  verv  well  with  tin 
measurements  from  experiment  lor  the  low  vele.  itv  impa.  t  .  as  , 
but  agreement  worsens  with  increasing  impact  velocity.  I'll  i  s  is 
to  be  expected  because  the  extent  of  fracture  damage  in.  reuses 
with  velocity,  and  fracture  was  not  accounted  for  in  these 
ca 1 cul at  ions .  Observed  crater  diameters  were  smaller  than 
computed,  possibly  because  the  occurrence  of  radial  cracks 
tends  to  relax  the  tangential  stresses  and  hence  reduce  the 
extent  of  plastic  flow.  The  measured  crater  depths,  on  the 
other  hand,  were  greater  than  computed,  suggesting  that  the 
occurrence  of  radial  cracks  facilitates  sphere  penetration. 

The  comparisons  in  Table  II  lend  credence  to  the  , ul.uluted 
stress  field  histories,  selected  samples  ol  which  are  displaced 
in  Figures  11  and  12. 

The  variation  in  tangential  stress  with  tinu  at  sevital 
distances  from  the  impact  crater  is  shown  in  Figure  11  tor  a 
60  m/s  impact.  Tin  tangential  stress  rises  s I ow 1 v  and  iu  aches 
a  maximum  at  about  8  ,s  after  impact,  and  then  de.  reuses  even 
more  slowly.  The  value  of  the  maximum  lulls  oft  stimuli  with 
distance.  At  large  distances  (three  crater  radii)  Iron  tin 
impact  axis,  the  tangential  stress  i.  initially  weak  I v  .  empress i vi 

Figure  12  indicates  the  variation  ol  peak  tangential  t  «  1 1  •  -i! 
stress  with  distance  from  t  lie  imp.nt  sit,  .  mitei  lot  three 


Table  11 


COMPARISON  OK  COMPUTED  AND  OBSERVED  RESULTS 
OF  PARTICLE  IMPACT  EXPERIMENTS 


IMPACT  VELOCITY 


5  m/s 

30  m/s 

_  60.  m/s 

Rebound  velocity  (m/s ) 

Observed 

4.8 

16.2 

30 

Computed 

4.9 

1  3.  1 

20.7 

Crater  diameter  (mm) 

Observed 

None  observed 

0.24  ±  0.01 

0.36  ±  0.01 

Computed 

0 

0.23 

0.4  7 

Crater  depth  (mm) 

Observed 

None  observed 

0.014  1  0.001 

0.034  ±  0.001 

Computed 

0 

0.012 

0.026 

31 


TANGENTIAL  STRESS  r„  (kbar) 


0  4  8  12 

TIME  (ps) 

V  A  1  1 

FIGURE  11  VARIATION  OF  TANGENTIAL  TENSILE  STRESS  V\ITH  TIME 
AT  1,  1  5.  2,  AND  3  CRATE  R  RADII  ui  FROM  IMPACT  SITE 


16.0 
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FIGURE  12  VARIATION  OF  PEAK  TANGENTIAL  TENSILE  STRESS 
WITH  DISTANCE  FROM  IMPACT  SITE  FOR  THREE 
IMPACT  VELOCITIES 


velocities.  Stress  rises  more  sharply  witli  distance  for  higher 
impact  velocities  and  falls  off  more  gradually.  Also,  the  peak 
stress  values  increase  with  impact  velocity,  even  though  the 
resolution  of  the  computer  printout  was  not  fine  enough  to 
capture  the  maxima  of  the  stress-distance  curves  in  Figure  12. 
Finally,  the  maxima  appear  to  occur  further  from  the  impact 
site  center  for  higher  velocity  impacts. 

The  computed  residual  stress  fields  produced  by  impact  at 
the  three  velocities  are  shown  in  Figure  13.  Compressive 
residual  stresses  are  predicted  near  the  specimen  surface  and 
extend  several  sphere  radii  from  the  impact  site.  A  residual 
tensile  stress  field  is  produced  beneath  the  impact  site.  The 
spatial  extent  of  both  fields  increases  with  impact  velocity. 
Maximum  values  of  residual  tensile  stresses  for  the  60  m/s  impact 
are  about  40  kbar  and  occur  approximately  one  half  sphere  radius 
beneath  the  impact  site.  Compressive  residual  stresses  reach 
maxima  of  about  2  kbar  for  the  30  m/s  and  60  m/s  impact  and 
occur  close  to  the  target  surface  at  about  one  sphere  radius 
from  tht'  impact  s  i  t  e  center. 

We  conclude  from  the  agreement  between  computed  and 
observed  crater  dimensions  that  the  stress  field  histories 
produced  by  hard  parti.  !e  impact  are  computed  reliably  and  in 
sufficient  letai  !  so  that  the  TROT  I  code  is  a  suitable  basis 
for  tin.  predictive  damage  capability.  In  Section  4.3  we 
construct  a  comput.it  ioiiui  mod.  I  tor  radial  cracking  around  impact 
sites  and  write  it  in  a  I  oi  suitable  for  incorpor.it  ion  into 
the  I'ROTT  code . 

*  •  1  foils  t  r  cic  t  i  on  ot  a  <'omput.it  ionul  1  r.u  lure  Model 

For  making  quantitative  predictions  ol  the  extent  of 
traeture  produced  in  a  target  impacted  hv  a  particle,  mat hema t  i  ca 1 
expressions  are  required  that  relate  true  ture  damage  to  the  causative 
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FIGURE  13  RESIDUAL  STRESS  FIELDS  IN  ZnS  AFTER  IMPACT  WITH  A  0  8-mm-riia 
TUNGSTEN  CARBIDE  SPHERE  AT  THREE  VELOCITIES 

Oj  ♦  Oy  +  n  -s 

- - in  kbar.  negative  values  ate  tensile 


Values  are  pressure 
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I  mat  In-mat  i  i  .i  I  express  i  ons  w. 


hav.  .h  v.  i  ■  ;  •  i  .  t  , ■  t  i  t  ..  -  a  i  .tapes  of  tin-  Ira-  tun-  prce.-ss 

a  Mil ,  !  iki  ii  :  h,  i  ,  -a  :  t  :  .<  !  i  a.  ■  t  u  re  mode 1  . 

I  '  i ■  i  r  >  :  i  i i '  ,  ,  :  ■ t  i  ,  tar.  m. nit  !  ".us t  ii« •  based  en  the  , i .  t  tia  1 
impac t -  induced  piivi,.,:  .  ,  n  t  ••  ...  .  urr  i  up.  in  tin  l.iru  t  that  lend  t . . 
tiu-  linai  i  t  j't  are  dam  n  >::;i  u,  - !  i  v  ■■!  tr.i  tan  in  ZnS  under 

single  parti'  1>  i  ••  pa.  :  wa  •  do. in-  e.i  i  .  t  i .  n  -  ■!  tia-  d.uapc 

prodn.  cd  h  v  imp. i  t  ,t  '.si  i .  m  ..  ! ,  .  it  i  <  ■  - ,  uid  is  d.  ..  rihed  in  the 
It'll  v'U  i  n subs,  lieu.  S  lbs,  ,  t  i . -u  •  •  ,1,  v,  'e;s  tin  expressions  tint 

describe  various  asp,,  t  -  ■!  t  in  p: u  n  m  <  n. 1 1  o p v  and  ilisiusses  the 

parameters  in  tins.  .  xpress  i.<a  ami  their  ral.  as  dynamic  Iraeture 
properties.  S.etien  ...  ,!  i  s.  uv-.s  i  n.  arpor.it  i  on  ot  tin  model  into 
tile  wave  propagation  .ode  d.s,  r  i  in  d  in  So.  tion  ..2,  and  recommends 
resoareh  tasks  lor  luture  programs. 

.  i.  1  Fhenotieiii)  1  ogv  ot'  Part  ic  le  Han.tve  in  ZnS 

Tin*  sequence  of  cracking  events  in  ZnS  was  determined  by 
direct  optical  observation  and  by  acoustic  emission  under  quasi- 
static  indentation  with  a  le.'  .spber.  '4  .aid  also  was  infer  cd  from 
post  test  examination  of  impacted  specimens.  ”  The  sequence  i  as 
follows:  As  the  sphere  is  pushed  into  tin-  ZnS  surface,  plastic 

flow  occurs  under  tile  contact  area  and  a  permanent  impression 
torms.  The  diameter  increases  with  increasing  load,  and  is  i  i  neat- 
in  a  lop-lop  plot.  Figure  10.  At  a  certain  threshold  load,  about 
five  roughly  oqui spaced ,  shallow  surface  cracks  form  and  spread 
radially  outward  from  the  crater  periphery  on  plains  perpendicular 
to  the  surface.  At  increasing  loads  these  cracks  increase  in 
length  and  slightly  in  depth,  maintaining  a  thumh-nai I -shaped 
profile  and  additional  radial  cracks  appear.  The  number  of  large 
(—200  am)  radial  cracks  attains  a  maximum  of  10  to  12,  whereas 
smaller  radials  increase  in  number  beyond  SO  and  maintain  an 

c  r .  1 1  ■  k  s  , 


exponential  size  distribution. 

At  loads  higher  than  those  required  to  form  radial 


oao  or  two  penny-shaped  cracks  containing  the  indentation  axis 
form  beneath  the  plastically  impressed  surface.  These  cracks, 
called  median  vents,  expand  radially  outward  and  often  join  with 
the  radial  c-acks. 

Still  higher  loads  are  required  to  produce  lateral  vents — 
cracks  that  form  and  propagate  parallel  to  the  indented  surface 
between  the  radial  cracks.  The  number  and  sixes  of  lateral  crack 
increase  with  load,  reaching  maximum  values  similar  to  those  of 
the  radial  cracks. 

As  the  load  is  removed,  all  three  tvpes  of  cracks  grow 
further.  The  depth,  and  to  a  lesser  extent  the  diameter,  of  the 
plastic  impression  recovers  somewhat,  but  residual  plastic 
strains  apparently  hold  the  radial  cracks  open.  Sectioning 
studies  show  a  heavily  deformed  or  finely  fractured  zone 
immediately  below  the  impression,  and  often  two  or  three  levels 
of  lateral  cracks. 

4.3.2  Rudiments  of  the  Fracture  Model 

The  observations  described  in  the  previous  section  suggest 
that  the  fracture  damage  produced  in  ZnS  by  particle  impact  occur 
by  the  nucleation  and  growth  of  cracks  in  planes  and  directions 
dictated  by  the  stress  field  histories  in  the  vicinity  of  the 
impact  site.  Thus,  the  computational  fracture  model  should 
consist  of  an  expression  describing  crack  nucleation  for  comput¬ 
ing  the  number  of  cracks  produced  by  a  given  impact,  and  an 
expression  describing  crack  growth  for  computing  the  size  of  the 
nucleated  cracks.  In  addition  it  is  necessary  to  treat  the 
relaxation  of  the  impact-induced  stress  field  that  results  when 
cracks  develop. 

Crack  nucleation  is  the  result  of  activation  of  pre¬ 
existing  flaws  during  passage  of  the  stress  pulse.  Flaws  (or 
weak  sites)  such  as  intergranular  or  transgranular  cracks,  pores, 
grain  boundaries  or  triple  points,  and  machining  and  grinding 


damage  art'  inherent  in  materials,  and  these  flaws  may  become 
unstable  and  propagate  when  traversed  by  a  stress  pulse  of 
sufficient  amplitude  and  duration.  Linear  elastic  fracture 
mechanics  relates  the  critical  stress  amplitude  to  the  crack 
length  by  the  stress  intensity,  and  shows  that  the  larger  the 
preexisting  crack,  the  lower  the  continuum  stress  required  to 
activate  it.  For  the  situation  where  the  load  is  not  quasistatic, 
but  rather  is  applied  as  a  short-lived  stress  pulse,  the  duration 
of  the  stress  pulse  must  he  considered  in  crack  instability 
ana  1  vses .  20 

Considerations  of  stress  intensity  histories  experienced 
by  cracks  whose  lengths  are  comparable  to  the  distance  a  stress 
wave  can  propagate  during  the  pulse  life  suggest  that  the  attainment 
of  a  critical  level  of  stress  is  necessary,  but  not  sufficient, 
for  crack  instability  A  sufficient  condition  is  that  the  stress 
level  be  attained  Cor  some  minimum  time.  Experimental  results  on 
several  polymers  and  structural  metallic  alloys  support  this 
premise.  21 

Thus  the  expression  for  crack  nuclcation  must  relate  the 
number  N  and  sizes  R  of  inherent  flaws,  the  resistance  of  the 
material  to  crack  activation  (taken  as  a  fracture  toughness 
appropriate  for  particle  impact  loading  rates,  and  the 

minimum  time  for  instability  tmjn  to  the  amplitude  o,  duration 
T0,  and  perhaps  the  shape  <5  of  the  stress  pulse.  We  can  write 
the  general  form 

dN/dt  -  f(N,  R,  Kjd,  train.  Tt>,  f)  (4-1) 
for  the  nuclcation  function. 

In  our  previous  work  on  ZnS ' ^  we  showed  that  the  shape  of 
the  flaw  size  distribution  is  exponential,  having  the  form 

\  =  No  exp(-R/Rn)  (4-2) 

where  Ng  is  the  cumulative  number  of  flaws  per  unit  volume  witli 
radii  greater  than  R,  Nf,  is  the  total  number  ol  flaws,  and  Kn 
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is  a  size  parameter.  The  quasistat ic  fracture  toughness,  Kjc , 
has  been  determined  in  previous  experiments ^2 , 23  to  |1(, 

0.75  i  0.01  MN  m-^^,  and  little  rate  effect  is  expected  in  this 
semibrittle  material. 

The  approximate  functional  form  of  Eq .  (4-1)  will  be 
established  in  future  work  by  correlating  the  measured  crack 
size  distributions  from  the  experiments  in  Reference  8  with  the 
computed  stress  field  histories  described  in  Section  4.2. 

Several  expressions  for  crack  growth  have  been  postulated 
from  theoretical  considerations .24-26  The  crack  size  distributions 
predicted  by  these  expressions  are  illustrated  in  Reference  8 
and  compared  with  measured  crack  size  distributions  at  impact 
sites  in  ZnS.  Since  all  theoretical  expressions  result  in 
distributions  similar  to  that  observed,  we  chose  the  simplest 
expression  for  the  present  fracture  model,  namely, 

dR/dt  =  kC£  (4.3) 

where  dR/dt  is  the  crack  growth  rate,  C£  is  the  longitudinal  wave 
speed,  and  k  is  a  constant  between  zero  and  1.  In  previous  work 
on  dynamic  brittle  rock  fracture ,^7  the  theoretical  value  ^ 
k  =  0.38  produced  good  agreement  between  computed  and  observed 
damage . 

As  fracture  damage  develops,  the  stresses  in  the  target 
relax.  When  radial  cracks  form  and  grow,  the  tensile  stresses 
on  their  surfaces  fall  to  zero  and  this  unloading  information 
is  communicated  to  the  surrounding  material  by  a  wave  that 
propagates  away  from  the  newly  formed  crack.  Tensile  stresses 
in  the  plane  of  the  crack,  however,  may  remain  high,  so  that  the 
stress  relaxation  process  may  be  highly  anisotropic.  Thus  the 
occurrence  of  cracking  modifies  importantly  the  stress  field 
history  produced  by  an  impacting  sphere,  and  may  explain  the 
increasing  divergence  with  impact  velocity  (increasing  cracking) 
of  computed  and  observed  crater  dimensions  and  rebound  velocities 
(Table  11).  Clearly,  a  reliable  predictive  capability  requires 
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that  the  effect  of  developing  fracture  damage  on  the  state  of 
stress  be  continuously  accounted  for  during  a  computational 
simulation. 

The  basis  for  the  treatment  of  crack-induced  stress 
relaxation  is  the  concept  of  a  two-component  system:  solid 
material,  and  void  inside  the  open  cracks.  The  specific  volume 
of  the  system  changes  when  loaded;  this  change  is  due  partially 
to  strain  in  the  solid  and  partially  to  change  in  the  void 
volume.  The  void  volume  calculation  follows  the  analysis  of 
Sneddon-®  for  a  penny-shaped  crack  in  an  elastic  material. 

During  the  wave  propagation  calculation,  the  following  sequence 
of  events  is  envisaged  while  the  material  is  in  tension: 

(a)  Initial  tensile  loads  cause  only  elastic  volume  changes 
in  the  solid  until  a  threshold  stress  (equal  to  the 
dynamic  fracture  strength)  is  reached. 

(b)  When  the  stress  exceeds  the  fracture  threshold,  cracks 
begin  to  nucleate  and  grow,  and  the  void  volume  produced 
by  the  cracks  acts  to  decrease  the  volume  change  required 
of  the  solid.  Thus  the  tensile  stress  in  the  solid  is 
lower  than  it  would  be  for  undamaged  material  under  the 
same  volume  change;  lienee,  the  effective  modulus  of  the 
solid  has  decreased. 

(c)  Witli  increasing  stress  and  continuing  volume  change',  a 

point  is  reached  where-  the  void  volume  increase  just 
equals  tile  applied  volume  change.  Here  there  is  no 
change  in  solid  volume  and  hence  no  change  in  tensile 
stress:  The  stress-volume  path  has  reached  a  peak  and 

tlie  effective-  modulus  is  zero. 

(d)  Witli  further  volume  changes,  the-  increase  of  void  volume 
(by  growth  and  nucleation  of  cracks)  tenuis  to  excc-cd  the' 
applied  volume  change1.  Then  the  solid  volume  change  is 
negative  and  the  tcuisilc1  stress  is  decreasing.  The 
effective  modulus  is  nc'gative  during,  this  period. 
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(e)  If  compressive  volume  changes  occur  at  any  time,  there  wi  1 
usually  be  a  decrease  in  both  void  volume  and  solid  volume 
Then  the  effective  modulus  is  positive  as  in  (a). 

(f)  If  no  net  volume  change  occurs,  there  will  usually  be  an 
increase  in  void  volume  and  a  corresponding  decrease  in 
solid  volume:  The  effective  modulus  in  such  a  case  is 
infinite . 

The  preceding  sequence  of  events  is  treated  in  the  present 
formulation  of  the  behavior  cf  material  undergoing  fracture.  A 
derivation  of  the  effect  of  damage  on  the  material  stiffness  is 
given  below  to  indicate  the  basis  of  the  method.  Here  the  case 
is  considered  in  which  flaws  are  present  but  no  nucleation  or 
growth  of  cracks  occurs:  The  resulting  modulus  is  that  which  is 
appropriate  for  a  residual  strength  calculation.  The  fundamental 
relation  is  that  the  total  volume  change  is  the  sum  of  solid  and 
void  volume  changes: 

AV  =  AVS  i  AVv.  (4-4) 

The  void  volume  change  is  derived  from  the  analysis  of  Sneddon  28 
for  the  opening  of  a  penny-shaped  crack  in  an  elastic  medium 
under  uniform  tension.  The  half-opening  of  the  crack  faces  is 

6  =  A(.l_-  y2h:o 
ii  E 


where  c  is  the  half  crack  length 

0  Is  the  stress  applied  normal  to  the  crack  plane 
E  is  Young's  modulus 
v  is  Poisson's  ratio. 


The  void  within  the  crack  faces  is  an  ellipsoid  with  semiaxes 
c,  c  and  6,  so  the  volume  is  ^ 

Vic  -  43^-2  =  (4_5) 


To  determine  the  total  relative  crack  volume,  a  ium  is 


Vv  =V  NiVir 

i 
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made  over  all  cracks 


per  unit  volume. 


where  is  the  number  of  cracks  of  volume  V^t. 
Then  F.q.  (4-4)  becomes 


r.v  =  Ao  Vs  16(1  -  v2) A(j 
V  (K  +  4u)V  3E 


2 


(4-6) 


where  Ao  is  a  change  in  stress 

K,  ij  are  buLk  and  shear  moduli 
V  is  the  total  specific  volume. 

Dividing  Eq.  (4-6)  by  Ao,  we  derive  an  effective  modulus  M. 


1  =  Va/X--  + 

M  K  +  Ap 


16(1 


Z  Nici- 


(4-7) 


This  compliance  is  similar  to  that  which  would  be  obtained  for  a 
composite  made  of  solid  material  plus  void.  The  term  K  +  (4/3);. 
is  the  stiffness  of  the  solid  and  Vs/V  is  its  volume  fraction. 

The  last  term  on  the  right  in  Eq .  (4-7)  must  be  interpreted 
then  as  the  compliance  of  the  cracks.  This  expression  shows 
that  the  compliance  of  the  composite  increases  with  increasing 
void  fraction,  that  is,  increasing  number  and/or  increasing 
size  of  cracks.  Since  both  the  strain  in  the  solid  and  the  void 
volume  are  tensor  quantities,  the  stress  relaxation  is  appropriate¬ 
ly  anisotropic. 


4 . 4  Future  Work 

The  final  steps  in  the  development  of  the  predii  .  !ve  capability 
for  particle-impact-induced  fracture  damage  are  to  evaluate  the 
nucleation  function,  Fq .  (4.1);  insert  the  fracture  model  (Section  4.1) 
into  the  TROTT  code  (Section  4.2);  simulate  computationally  several 
selected  particle  impact  experiments;  compare  the  computed  and  observed 
fracture  damage ;  and,  if  necessary,  modify  the  expressions  and  para¬ 
meters  in  the  model  to  attain  good  agreement. 

The  nucleation  function,  Fq .  (4.1),  should  he  evaluated  by 
correlating  computed  stress-strain  field  histories  for  experiments 
of  several  impact  velocities  with  measured  radial  crack  sizi 
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distributions.  We  expect  these  correlations  to  be  difficult  because 
of  the  high  gradients  in  stress  and  strain  near  the  impact  site.  More 
suitable  experiments  are  flat  plate  impact  experiments  such  as  those 
reported  in  Reference  29,  because  the  specimens  contain  relatively 
large  volumes  of  material  that  experience  nearly  constant  stress  and 
strain  histories.  Thus,  statistical  data  on  the  number  and  sizes  of 
cracks  can  be  obtained  for  correlation  with  a  stress  or  strain  state. 
Limited  data  on  ZnS  exist  ^9  anc]  should  be  examined  to  help  in  the 
correlations. 

The  fracture  model  presented  in  Section  4.3  contains  only  one 
orientation  of  cracks  and  treats  only  an  exponential  crack  size 
distribution.  The  model  should  be  extended  to  pertain  to  cracks  in 
three  orthogonal  directions  so  that  radial,  ring,  and  lateral  cracks 
can  be  treated  simultaneously.  Furthermore,  since  nonexponential 
crack  size  distributions  have  been  observed  after  considerable  crack 
growth,  appropriate  mathematical  descriptions  of  the  distributions 
should  be  developed  and  incorporated  into  the  model. 

These  modif icaCions  have  been  studied  extensively  to  determine 
their  feasibility  and  possible  implementation  procedures.  The  use 
of  three  orientations  of  cracks  introduces  significant  convergence 
problems  in  the  stress  and  damage  calculat ional  algorithm.  Currently, 
with  one  crack  orientation,  an  iteration  procedure  is  used  to  determine 
one  unknown  strain  from  the  nucleation,  growth,  and  crack-opening 
functions.  With  three  orientations,  an  iteration  procedure  for 
three  unknowns  will  be  required.  Because  of  the  nonlinearity  of 
the  functions,  convergence  of  the  iterations  will  be  difficult  to 
obtain,  but  techniques  developed  at  SRT  International  for  similar 
problems  may  be  applicable. 

Representation  of  more  complex  crack  size  distributions  can  be 
done  in  several  ways.  One  method  we  have  used  in  other  cases  is  to 
describe  the  distribution  by  a  series  of  N-K  values  instead  of  by  an 
analytic  function.  This  method  is  straight f orward  to  use,  but  it 
adds  greatly  to  the  computer  storage  required.  An  alternative  method 
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is  to  use  a  higher  order  function  to  represent  the  ir.uk  size  ilist  r  ibu- 
tion;  this  method  requires  derivation  <>i  nueleation  and  growth  laws  t  <  >  r 
each  of  the  parameters  of  the  distribution.  A  third  alternative  i •  .1 
combination  of  tile  above:  Describe  the  crack  size  distribution  by  a 
series  of  exponential  functions,  each  governing  a  narrow  range  of 
crack  radii.  These  three  possibilities  should  be  explored  tor  their 
appropriateness  to  the  problem  and  the  efficiency  with  which  thev  can 
be  executed. 

These  modifications  should  be  done  concurrently  with  work  to 
correlate  the  nueleation  and  growth  expressions  with  observed  damage . 

4 . 4  Summary 

Considerable  progress  has  been  made  in  developing  a  predictive 
capability  for  fracture  damage  in  ceramics  produced  by  particle  impact. 

A  two-dimensional  l.agrangian,  finite  difference,  wave  propagation 
code,  TROTT,  was  provided  with  a  slide  line  capability  and  used  to 
simulate  the  impact  of  an  800-i.m-dia.  WC  sphere  and  a  ZnS  plate  at 
velocities  of  5,  10,  and  60  nt/s.  Plastic  flow  was  allowed  to  occur 

according  to  an  experimentally  determined  stress-strain  relation,  but 
cracking  and  the  associated  stress  relaxation  were  not  treated  in 
these  initial  calculations.  The  computation  simulated  qualitatively 
all  of  the  important  aspects  of  a  particle  impact  event,  and  gave 
good  quantitative  agreement  with  experimental  measurements  of  rebound 
velocity  and  crater  dimensions. 

Mathematical  expressions  that  describe  crack  development  in  X.nS 
under  particle  impact  load  histories  wore  deduced  from  observations  of 
the  number  and  sizes  of  radial  cracks  produced  by  impact  at  various 
velocities.  These  expressions  treat  damage  development  as  a  crack 
nueleation  and  growtli  phenomenon  and  constitute  .1  computational  fracture 
model  that,  written  as  a  subroutine,  can  he  used  in  conjunction  with 
the  wave  propagation  code  to  calculate  the  extent  ot  tract  lire  damage 
resulting  from  a  given  particle  impact.  The  parameters  in  the  model, 
because  they  are  materia l-spec if ic  and  In  in  <  act  as  dynamic  fracture 
properties,  can  be  used  to  rank  and  evaluate  i.mdidati  materials  tor 
various  impact  erosion  app  1  i  1  a  t  i  oils . 
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Objectives  i>t  future  work  arc  to  establish  the  I  inn  t  i 
he  eraek  iniclcation  I'xpression,  tu  incorporate  tin  tr.ietni 
he  wave  propagation  rode,  and  to  siniulate  e  ompn  t  a  L  i  ona  1  1  v 
mi t i *  le  impart  experiments. 
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